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The real part of the self-energy of interacting two-dimensional electrons has been calculated 
in the i-matrix approximation. It is shown that the forward scattering results in an anomalous 
term leading to the vanishing renormalization factor of the one-particle Green function, which is 



a non-perturbative effect of the interaction U. The present result is a microscopic demonstration 
of the claim by Anderson based on the conventional many-body theory. The effect of the 
damping of the interacting electrons, which has been ignored in reaching above conclusion, has 



£H , been briefly discussed 

o 
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The nature of the low-energy excitation of the two-dimensional system is of great interest re- 
cently. The Fermi-liquid picture was considered to be valid from diagramatic studies ,EP while it 
was suggested by Anderson that the anomalous behavior of the forward scattering phase shift leads 
to the renormalization factor Z = 0, i.e., the breakdown of the Fermi- liquidBu) This remarkable 
suggestion has attracted much interest,a , H' and several calculations of the self-energy of the two- 
dimensional Hubbard model have been carried out based on the ^-matrix approximation jQ'Q'&tJ) 
T3 ; 

in which the self-energy is approximated by the summation of ladder diagrams of the particle- 

O ■ 

particle process. In these calculations, however, only the imaginary part of the self-energy has been 

considered, and the real part has not been studied in detail. 

In this paper, we will calculate explicitly the real part of the self-energy of the two-dimensional 

h : 

Hubbard model by the i-matrix approximation, and show that the singularity of the i-matrix in 
the forward scattering region gives rise to an anomalous term to the real part of the self-energy, 
which leads to the renormalization factor Z = 0. This result is in accordance with the claim by 
Anderson. 
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We consider the asymptotic behavior of the self-energy S(fc, e + i<5) in the limit of \k — k-p\ ^ 
^F; H *C £Fi &f and £f being the Fermi momentum and the Fermi energy. In general, the shape 
of the Fermi surface near the point k in the momentum space can be approximated as a parabolic 
curve, as is schematically shown in Fig. 1(a) and 1(b). Linearizing the energy dispersion in the 
normal direction of the Fermi surface, we assume the following energy dispersion.!!!^ 

Ck = e k - H = v k x + -^-ky, (1) 

where we take the origin of momentum at the point nearest to k on the Fermi surface, as is 
shown in Fig. 1(b), /x, vq and Ag are the chemical potential, the Fermi- velocity and a constant, 
respectively. Here, the momentum dependence of the velocity is neglected. We calculate the 



Fig. 1. The local shape of (a) the Fermi surface can be approximated by (b) a parabolic curve, (c) The energy e 
dependence of the self-energy in the vicinity of vk is focused on. 



self-energy £(fc, e + i<5) at the point k = (k, 0). In the Hubbard model, the ^-matrix is given by 

n«'+w= T*miTm (2) 

*(,,* + *) = £-: -^±22 -,, (3) 

k Zk+q/2 + C,-k+q/2 ~ X — 10 

where K(q, x + iS) is the particle-particle correlation function. Noting that K(q, 0) as a 
function of q is regular in the forward scattering region, q ~ 2kp, we introduce K$ = 
linig_ > 2fc P [h m z^o K(q, x + i5)], which has a contribution from the high-energy region, and reflects 
the whole band structure. In contrast to this, q- and x-dependences of K(q, x + i5) — Kq in the 
region of \q — 2kp\ <C /sp, \x\ <C eF reflect details of a scattering process near the Fermi-energy. 
Hence, we write the i-matrix as, 

T (*,x + -<6)= 1 + UielK - q U * +iS) _ Ko] (4) 

u " s TTuK a (5) 



We will focus only on the contribution of the t-matrix to £(fc, e + i<5) from the forward scattering 
region, i.e., the total momentum of the particle-particle ladder is nearly 2k. In the following 
calculations, we introduce momentum cut-offs k c and k' c for x and y components of the momentum 
in the energy dispersion, and momentum integrations for intermediate states of a scattering process 
are carried out within these cut-offs. 

In the i-matrix approximation, the expressions for the real and imaginary parts of the self-energy 
are given by 

f d 2 q f°° dx 
ReS(fc,e + id) = / - — r^ / — ReT(q, x + id)lmGo(q — k,x — e — i<5) 

J {2tT) z Je VT 

/d 2 <7 f°° dx 
7^2 I —lmT(q,x + i5)ReG (q-k,x-e-id) (6) 

ImE(fc, e + iS) = - f -^ f —JmT(q, x + id)ImG (q -k,x-e-id). (7) 

J (2ir) z Jo it 

We will evaluate the asymptotic forms of Re£(fc, e + id) and ImS(fc, e + id) under the condition of 
|e|, \vk\ <C e c , where e c is defined as Min(vok c ,Aok' 2 ). Dividing x integration in [e,oo] in the first 
term of eq. (6) into [0, e] and [0, oo], we introduce Re£i(fc, e+id), ReE 2 (fc, e+id) and Re£3(fc, e + id) 
as 

Re£(fc, e + id) = ReEi(fc, e + id) + Re£ 2 (fc, e + id)+ ReS 3 (fe, e + id) (8) 

/d 2 o r°° dx 
—-^ / — [ReTfg, x + id)ImG (q -k,x-e-id) 
(2lty Jo 7T 

-lmT(q, x + id)ReG (q -k,x-e- id)] (9) 

/d 2 a f e dx 
— 4r / — ReT(q,0)lmG (q-k,x-e-id) (10) 

(2tt) z Jo vt 

f d 2 o f e dx 
Re£ 3 (fc,e + W) = - / — -^ / — [ReT(q,x + id) - ReT(q,0)l ImG (q - k,x - e - id). (11) 
J (2tt) z Jo vt 

It is seen that ReEi(fe, e + id) is a function of only e — vk and ReX 2 (fc, e + id) a function of only e 
for the dispersion given by eq. (|]). Reflecting the regular behavior of ReT(q, 0), Re£ 2 (fc, e + id) is 
a regular function of e. Actually for |e|, \vk\ <C e c , we obtain 

Re£i(fc,e + i<5) ~ ReS(0,0) + ci(e - uJfc) (12) 

Re£ 2 (fc,e + i5)~c 2 e, (13) 

where ci and c 2 are constants. On the other hand, we note that ReS3(fc, e + id) and ImE(fe, e + i<5) 
are real and imaginary parts of a function which has a contribution from only the low-energy region, 
i.e., \q x \ < fe c , |^| <C fc£, |x| < e c . 

ReS 3 (fc,e + i(5)+iImS(/c,e + W) (14) 

-- t T 1 / V 7T f — [T(q,x + i5)+^ ff ]ImGo(q-fc,x-6-i<5), (15) 

J-A <57T J_A' *7T JO 71" 



where the cut-offs A and A' are small values compared to k c and k' c , respectively, and U c s = 
— lim q ^o[lim x ^oT(q,x + i5)]. We introduce S s (fc, e + id) by subtracting the on-shell value of 
Re£ 3 (fc,e + i(T), 



E s (fc,e + i<5) = [ReS 3 (fc,e + i<y)-ReE 3 (fe,uA; + W)] + iImE(fc, e + i<5). 



(16) 



ReS s (fc,e + id) is of interest to us, while ReE(fc,e + id) — ReS s (/c,e + id) = ReEi(fc,e + id) + 
ReS2(fc,e + i<5) + HeT,^(k,vk + i<5) is related to various renormalizations such as the shift of the 
chemical potential, the renormalization of the Fermi velocity, the effects of which are taken into 
account by replacing vq and Ao in eq. (|l]) and the renormalization factor Zq of the Green function 
Go(k,e + i<5). Then, we obtain the renormalized self-energy T,*(k,e + i<5) and the renormalized 
Green function G*(k, e + i<5) as 



G*{k,e + id)~ Z " 



(17) 
(18) 



(e-?;A;)-£*(fc,e + iJ) 

where G$(k, e + i<5) = Z^jie — £|) and £| = u /c^ + A^/2. The subscript Go — ► Gq on the right hand 
side (r.h.s.) of eq. (17) indicates that the calculation of E s (fc,e + i<!>) should be performed using 
the renormalized Green function, Crg(fe, e + id). We note that the on-shell value of ReE3(fc, e + iS), 
KeT,^(k,vk + id), is of the order of (vk) 2 , and that this term gives the renormalization of a quasi- 
particle energy only of the order of k 2 . 

Next we need to estimate S*(fc,e + id). Under the conditions of \q x \ <C k c , \q y \ <C k' c , \x\ <C e c , 
the asymptotic form of K*(q, x + id) — K$ is given by 



iZlx 



K*(q,x + id)-K* 



4ttA 1 / 2 v(x - vq x - Aq 2 /^) 1 / 2 

Zq x 

A7rA 1 / 2 v(-x + vq x + Aq 2 /Ay/ 2 



(x > vq x + Aq 2 y j\) 
(x<vq x + Aq 2 /4), 



(19) 



where the particle-particle correlation function K*(q,x + id) is calculated using the renormalized 
Green function G^k, e + id). Using these expressions, we obtain the asymptotic form of ReS*(fe, e + 
id) in the limit of [e — vk\ <C A~ x \J 2 ftZ§k 2 for the case of e > vk, for example, as 



Re£*(fc,e + i<5) 



4vr 2 



V J-X 



dq % 



dxRe 



-u cS 



-u cS 



1 + 



iU cS Z Q x 



4irA 1/2 v(Aql/4 + e - vk) 



T72 



1 + 



iU cS Z x 



^A l/2 v{Aq 2 y /4) 1 



~ 



(e — vk) 



irZt 



Jo 



dq y 



q y 



2 \ 1/2 

+ 1 arctan 



U c sZ e 



47rA 1 /2 u ( e _ w jfc)i/2 ^2/4 + i 



1/2 



v\ 1/2 

-r- arctan 



U c sZ e 



K 4nAV*v{€ - vk) 1 / 2 (gJ/4) 



1/2 



sgn[k]- 



vk) 



~ sgn[fc] - 



vk) 



dq y 



log v, 



% 



1 \ 1/2 

+ 1 



2\ V 2 ' 



where 



U cS Z 2 \k\ 



>1. 



(20) 



(21) 



A 1 / 2 ^ - vk) 1 / 2 
The calculations for the case of e < vk are similar and we obtain the final result in the limit of 

|e-ufc| <C A^U^Zfe 2 as 



-r, «/, .px (e-vk), A\e — vk\ 

ReT.*(k,e + i5) ~ c- — - — - log 



Za 



U^k 2 ' 



(22) 



where 



- (e > vk,k > 0) 




1 

2 
1 



(23) 



(e <vk,k> 0) 

(e > vk,k < 0) 

(e <vk,k < 0) 

Reflecting the existence of a logarithmic singularity, the renormalization factor Z = lim e _>„A:(l — 
<9ReX/9e) -1 exhibits an anomalous behavior, Z = 0, except in the limit of e — > vk — in the case 
of vk > 0. 

In addition to this, we have evaluated the asymptotic form of Im£*(/c, e + i<5) in the same limit 
of |e — vk\ <C A^^^U^Z^k 2 , which is given by 

ImS*(fe, e + i<5) ~ Im£*(fc, vk + i<5) 



+ 



(e — vk) 
Zo 



1_ Ae-vk\ 



+ a log 



.A|e — f k\ 



+ b 



t wt ua.^ U ctt Z ok 2 1 ^ cff Z 2 |fc| 

ImS (k,vk + lo) ~ ^ — log — _ ,„ " V = Ifc, 



(24) 



(25) 



^ ~° AV 2 el' 2 

where a and b are constants which are different in each region of e > vk and e < vk. Especially 
for e — ► wfc, only the first term on the r.h.s. of eq. (24) survives. So far, the real and imaginary 
parts of the self-energy have been examined, but E(fc, e + i<5) must be an analytic function of 
e. Actually, £*(£:, z) is given as follows for a complex variable z (Im z > 0) in the limit of 
\z - vk\ < A-if/^Z^ 2 . 

E*(jfe,*) =silm£*(fc,vife + i<$) 



(z — vk) 

Z~n 



i / A(z - vk) 

^[ log li^zW, 



, A{z-vk) \ 

+a V og ^zw) +h 



(26) 



where a' and b' are complex constants related to a and b in eq. (24). From the analyticity 
condition, T,*(k, e + id) is interrelated from e > vk to e < v k through the upper-half complex plane 
of e (Fig. 1(c)), and therefore the constants a' and b' in each region are mutually related by a simple 
relation. When we see £*(£;, z) as one analytic function of z, the logarithmic term in the real part 
is related to the second term of the r.h.s. of eq. (24) in the imaginary part. Actually, we obtain 
Re a' = —1/2 (1/2) for vk > (vk < 0). The leading term of the imaginary part when e ~ vk, 
however, is not the second term in eq. (24) but the first one. On the other hand, in the case of 
k = kp, ImS(A;F, e + i<5) ex e 2 log e was obtained J1P which corresponds to the first term on the r.h.s. 
of eq. (24), and was the basis of the claim of the Fermi- liquid state. As has been demonstrated, 
however, the present logarithmic term of the real part cannot be deduced from this imaginary part 
using the Kramers-Kronig transformation. This has not been noted in previous studies. 

Our new result of a logarithmic singular behavior in the real part of the self-energy implies that 
the renormalization factor Z vanishes in the low-energy limit. 

In addition, we note the relationship between the results in the i-matrix approximation and the 
second-order perturbation theory. When we focus only on the forward scattering process, the self- 
energy as an analytic function of complex variable z is obtained by the second-order perturbation 
theory aaL2P 

E*(fc,s)ociCfVlog ( *~' ,fc) , (27) 

and there is no logarithmic term in the real part of the self-energy In the case of the t-matrix 
approximation, we found that the self-energy is a logarithmic function not only of e — vk but also of 
U e g, as shown in eq. (26). Hence, the logarithmic term in the real part of the self-energy obtained 
by the i-matrix approximation is a non-perturbative effect of U, which implies that U = can be 
a singular point in two dimensions. 

Investigating the derivation in eq. (20) of the logarithmic term in the real part of the self-energy 
near the on-shell region, we can see that the logarithmic term is related to the property of the 
t-matrix in the region of \x — vq x \ < A^^^U^Z^k 2 , Aq 2 < A^ x \J 2 ^Z^k 2 , < \x\ < |e|. This indicates 
that the on-shell electrons are responsible to the logarithmic term of the self-energy. 

Anderson proposed that the special feature of the forward scattering process leads to the van- 
ishing renormalization factor, Z = 0, i.e., the breakdown of the Fermi-liquid, and our results 
have confirmed this based on the conventional many-body perturbation theory using the Feynman 
diagrams. 

So far the calculations of the self-energy have been carried out based on the Green function 
without damping, T^- The resulting logarithmic behavior in the real part of the self-energy, however, 
has been seen only in the energy region |e — vk\ <C A~ 1 U 2 s ZQk 2 Ti : , implying the importance of the 
self-consistent treatment. Actually, Tfc results in finite Z even when X — > 0, i.e., G(k,e + i<5) — 
Zo/[a(X)X + Hfc] with a(0) ^ 0. However, a(X) will have a strong dependence on X = e — vk as 



a(X) = ao[l + 0(|X|/r^)], which is different from a conventional Fermi-liquid. 
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